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Abstract 

Ultra-high energy cosmic rays are mostly charged particles and they are 
therefore deflected by magnetic fields on their path from their sources to 
Earth. An interesting phenomenon arising from these deflections is the ap- 
pearance of multiple images of a source, i.e. cosmic rays with the same 
energy coming from the same source that can arrive to the Earth from dif- 
ferent directions. In this work we present a technique to identify secondary 
images, produced by the regular component of the galactic magnetic field, 
benefiting from the fact that near caustics the flux is significantly magnified. 



1 Introduction 

Galactic and extragalactic magnetic fields do not only deflect the trajectories of 
charged cosmic rays but they can also amplify or demagnify the flux arriving from 
a source, modifying its spectrum. A detailed analysis of magnetic lensing effects 
can be found in refs. |[T]|2l[3]ffl- An interesting phenomenon caused by magnetic 
lensing is the existence of multiple images, i.e. cosmic rays of the same energy 
coming from the same source that can arrive to the Earth from several different di- 
rections. The secondary images appear in pairs in the same direction of the sky at 
some critical energy Eq which depends on the magnetic field model, the cosmic ray 
charge Z and the direction to the source. When a pair of new images appear in the 
sky their flux is strongly magnified. For energies lower than the critical energy Eq, 
the arrival directions of cosmic rays of the secondary images move away from each 
other and their magnification decreases. Typically the flux of particles arriving to 
the Earth from the two secondary images in the energy range [0.9Eq, Eq] is an 
order of magnitude larger than the flux that would arrive if no lensing effects were 
present. Therefore a clustering of events with similar energy is expected around 
the direction where the secondary images appear and it is thus possible to benefit 
from this fact and design a method to search for such signals in cosmic ray data as 
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proposed in this paper. 

In analogy with the gravitational lensing effect, we call caustics the positions of 
the sources for which a pair of images with diverging amplification appear. These 
form closed lines in the sky whose location is a function of the energy. The cor- 
responding positions of the pairs of new images are the critical lines. Their shape 
and position as a function of the energy are characteristic of each galactic mag- 
netic field model. The identification of the direction and energy at which the new 
images appear would then provide new insights on the galactic magnetic field. An- 
other handle on magnetic fields can be obtained from the observation of ultra-high 
energy cosmic ray multiplets (events with different energy that come from a single 
point-like source), which could be used to determine the integral of the perpendic- 
ular component of the galactic magnetic field along the cosmic ray trajectory [|5]. 

We note that for the galactic magnetic field model we consider the first caustics 
appear at energies Eq/Z ps 30 EeV. If the cosmic ray composition was heavy, the 
lensing effects we study will only appear at energies above the maximum energies 
detected by existing experiments. Hence, our results essentially apply to proton 
cosmic ray sources. In addition, even in the case in which the source produces 
heavy nuclei with a unique value of Z, photodisintegration during propagation will 
lead to a distribution of values of the charge reaching Earth, and the caustics will 
appear at different energies for each charge, diluting the possible signals. Effects of 
magnetic deflections of heavy nuclei in galactic magnetic fields have been recently 
discussed in 

The paper is structured as follows: in section 2 we describe the general known 
characteristics of the galactic magnetic field and review some magnetic lensing 
facts. In section 3, we explain the method used to search for signals of multiple 
images of a source and to characterize them. In section 4, we apply the method 
to simulations and estimate the density of sources for which the effect is expected 
to be observable. Finally, section 5 is a summary of the results and contains our 
conclusions. 



2 Galactic Magnetic Field and Magnetic Lensing 

Although there are several observational methods to obtain information about the 
galactic magnetic field (Zeeman splitting, polarized thermal emission from dust in 
clouds, polarization of starlight, synchrotron radio emission and Faraday rotation 
of polarized sources (7]|8]|9]), this field is not yet precisely known. 

The galactic magnetic field has a large-scale regular component and also a tur- 
bulent one. From the analysis of the polarization of starlight it is known that the 
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local regular component is approximately parallel to the galactic plane and fol- 
lows the local spiral arms iTTOl . The regular component has a local value between 
B reg ~ 2-4 fiG according to Faraday rotation measures iflXI or synchrotron mea- 
surements [8] (the possible origin of this difference is analyzed in Ifl2l [TBI [T41 ). 
Reversals in the direction of the magnetic field between neighboring arms have 
been detected with Faraday rotation data, which provide evidence that in the inner 
arm magnetic fields are counterclockwise when viewed from the North Galactic 
pole, while in the local region the fields in the interarm regions are clockwise Q. 
The galactic halo field remains a puzzling subject since rotation measures indicate 
that it appears to be antisymmetric with respect to the galactic plane [7 ], but it is not 
yet understood how can a halo field with odd parity be superimposed with a disk 
field with even parity (H. The turbulent component has a root mean square ampli- 
tude of B rrns ~ (l-2)B reg and a typical coherence length L c ~100 pc lfT31 [Toll . 
The dominant effect on the deflection of high energy charged particles traveling 
through the Galaxy is produced by the regular component because it is coherent on 
scales much larger than L c . 

Even though the presence of extragalactic magnetic fields affects the propaga- 
tion of cosmic rays, they are probably not the predominant effect in the deflection 
of ultra-high energy charged particles. For instance, from the results in ifTTl one 
finds that for 10 EeV protons the deflections are smaller than 1 degree in 98% 
(90%) of the sky for distances to the source of 35 Mpc (70 Mpc). However, the 
strength of the extragalactic magnetic fields is quite uncertain and for example in 
|[T8l larger values for the deflections are found. We will only consider hereafter the 
deflection produced by the galactic magnetic field. 

In this paper we consider an illustrative galactic field model that includes some 
of the observed features, such as field reversals, and reproduces the local field 
strength. The regular component of the galactic magnetic field is modeled with a 
bisymmetric field with even symmetry (BSS-S) with structure and strength very 
similar to those used in [19] but smoothed out as described in (JJ. In this model 
the galactic magnetic field reverses its sign between the arms of the Galaxy and the 
field is symmetric with respect to the Galaxy's mid-plane. The local value of the 
field is taken as 2 /xG. For the dependence on z (i.e. perpendicular to the galactic 
plane) a contribution coming from the galactic disk and another one from the halo 
are considered: 

B reg (x, y, z) = B reg (x, y,z = 0) (- 1 + 1 ) , (1) 

\2cosh(z/,zi) 2 cosh(z / z 2 ) J 

with Z\= 0.3 kpc and z 2 = 4 kpc. 

The deflections caused by the magnetic field can lead to lensing phenomena, 
such as the energy dependent magnification and demagnification of the flux that 
modifies the energy spectrum of the source, or to the appearance of secondary 
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BSS-S E/Z=20 EeV 
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Figure 1: 'Sky sheets': directions of incoming cosmic rays in the halo that corre- 
spond to a regular grid of arrival directions at Earth, adopting the BSS-S magnetic 
field configuration, for particles with E/Z = 20 EeV (a) and 10 EeV (b). 



images EIEIO. The magnitude of the flux amplification depends on the arrival di- 
rection, on the ratio between the energy and charge E/Z of the cosmic ray and on 
the magnetic field model. The formation of multiple images and the flux magnifi- 
cation can be understood pictorially plotting for a regular grid of arrival directions 
at Earth the corresponding directions from which the particles arrived to the galac- 
tic halo. This is shown in Figure Q] for the model of the regular magnetic field 
configuration considered and for particles with E/Z = 20 EeV (top panel) and 10 
EeV (bottom panel), where 1 EeV = 10 18 eV. One may picture this distorted im- 
age of the sky seen from the Earth as a sheet (the 'sky sheet') that can be stretched 
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(a) (b) 

Figure 2: Position of the images in galactic coordinates for a source located at 
(1, 6)=(4.7°, — 27.4°)(a) and their amplifications as a function of the energy (b). 
Black symbols correspond to the main image while blue and red symbols corre- 
spond to secondary images. A star indicates the position of the source. Plus signs 
correspond to events with energies between 6 and 9 EeV, asterisks correspond to 
energies between 9 and 11 EeV, circles correspond to energies between 11 and 13 
EeV and squares correspond to energies above 13 EeV. 



and folded. A source located within a fold will have multiple images, i.e. cosmic 
rays of the same energy can arrive from this source to the Earth along several dif- 
ferent directions. At each energy the position of the folds in the plots indicate the 
location of caustics. If a caustic crosses the position of a source at an energy Eq, 
a pair of new images appears and their magnification diverges as 1/ \/Eq — E for 
energies close to Eq [2]. Therefore, a clustering of events with energies close to Eq 
is expected, increasing the probability of detecting the multiple images. As the flux 
magnification depends on the observed direction, at lower energies each secondary 
image will have different amplification. Pictorially speaking, the flux coming from 
a source in a region where the sheet is stretched will appear demagnified while 
that from a source in a compressed region will appear magnified. As the sky sheet 
changes with the energy, the spectrum observed from a given source is different 
from the emitted one HI. In Figure |2] we show for illustration the position of the 
images and their amplifications as a function of the energy for a source located at 
galactic coordinates (/, 6)=(4.7°, -27.4°). In this example, one has E /Z = 13.0 
EeV. 

Magnetic lensing phenomena also appear for turbulent fields 0. In this case, 
multiple images appear below a critical energy E c such that the typical transverse 
displacements among different paths after traveling a distance L in the turbulent 
field become of the order of the correlation length of the random magnetic field 
{firms ~ L c /L). This critical energy is typically higher than the corresponding 
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one for the regular magnetic field, but the folds produced by the random field near 
the critical energy are on a much smaller angular scale. For decreasing energies, 
the fraction of the sky covered with folds increases, however the magnification 
peaks become increasingly narrower and for E < E c /3 their integrated effect be- 
comes less noticeable. This can be pictured in terms of the 'sky sheet' (Figure [3]>, 
where small folds due to the turbulent component appear on top of the large folds 
produced by the regular component, being hence these last ones those giving the 
dominant effects. In Figure |4j we plot the same example as in Figure |2j but with the 
events propagated through both regular and turbulent components of the galactic 
magnetic field. We see that the main properties of the magnification peaks which 
are relevant for our analysis (their location, Eq and their approximate shape within 
10% of Eq, as well as the elongated shape of the images near the caustic) remain 
essentially unchanged. Hence, on the following we will focus on the effects of the 
regular field alone. 



E/Z=10 EeV 
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Figure 3: 'Sky sheet' adopting the BSS-S magnetic field configuration and in- 
cluding a turbulent component, for particles with E/Z = 10 EeV. The turbulent 
component was simulated with a Gaussian random field with zero mean and root 
mean square amplitude of B rms = 2/iG and a coherence length of L c = 100 pc. 

Let us note that a possible background to the signals due to the appearance of 
multiple images (i.e. a clustering of events with similar energies) could be pro- 
duced in bursting models for the sources of ultra-high energy cosmic rays. Due to 
the energy dependent delay in the arrival time of cosmic rays in the presence of in- 
tergalactic magnetic fields, this would imply that at a given time only cosmic rays 
with the same time delay, and therefore same energy, would be observed at Earth 
coming from a bursting source [20 1. However, the lensing features should appear 
only below some critical energy (at which the first caustics appear in the sky) while 
those due to bursting sources should be relevant up to the highest energies. 
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Figure 4: Position of the images in galactic coordinates for a source located at 
(1, 6)=(4.7°, —27.4°) (same as the example in Figure 2) (a) and their amplifica- 
tions as a function of the energy (b). Full symbols are used for events propagated 
through both the regular and turbulent component, while empty symbols are used 
to indicate the events which have only been propagated through the regular field. 



3 Method 

The purpose of this work is to design a method to search for signals of multiple 
images of a source, taking into account that cosmic rays from that source should 
appear clustered in the sky at energies close to the critical one. In order to de- 
fine what would be the typical signal of the clustering due to the appearance of 
multiple images of an extragalactic point-like source, we simulate sources that are 
randomly distributed in the sky and keep those that are crossed by caustics for ener- 
gies above a given threshold. We note that for the BSS-S magnetic field model that 
we adopted, sources located in approximately 24% of the sky would have caustics 
at energies above 10 EeV and assuming a proton composition, while above 20 EeV 



this fraction is only 7%. In Figure 5(a) we show the location of the one hundred 



randomly located sources from which we simulate the events which are shown in 
Figure |5(b)| 



For each source we simulate cosmic rays which we then propagate through 
the galactic magnetic field. These simulations are done backtracking antiprotons 
leaving the Earth up to a distance where the effect of the galactic magnetic field be- 
comes negligible fl}. The initial direction of the backtracked antiproton (that will 
be associated to the arrival direction of protons of the same energy) is recursively 
adjusted until the final direction points to the source with an accuracy better than 
10 -4 degrees. 
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Figure 5: One hundred simulated sources (asterisks) that have multiple images in a 
projection of the celestial sphere in galactic coordinates overlaid with a 'sky sheet' 
of the BSS-S galactic magnetic field model for particles with E/Z = 10 EeV (a) 
and ten events from secondary images from each of these sources (b). The radius 
of the circles is proportional to the energy of the events. 



We adopt an E~ 2 power spectrum at the source which is then modified when 
observed at the Earth by the magnification of the flux as a function of energy in the 
corresponding direction of the sky. We display in Figure |5(b")1 ten events from the 
secondary images of each source, with energies between 0.9 Eq and Eq. The en- 
ergy interval is chosen so as to take into account a reasonable experimental energy 
resolution, which is larger than the width of the amplification peak (defined as the 
energy above which the amplification is larger than 3) that is equal or smaller than 



8 




Figure 6: Distribution of events within a certain 'width' and 'length' in the tangent 
plane of the secondary images with energies between [0.9 Eq, Eq] for the one 
hundred simulated sources, with (black solid line) and without (blue dashed line) 
angular resolution effects (of 1 degree standard deviation). 



0.02 Eq for 50% of the simulations. Note that being these magnification peaks so 
narrow in comparison with the typical experimental energy resolution, it is unfea- 
sible to study detailed correlation between the location and the energy of the events 
appearing in these magnification peaks. The arrival directions of cosmic rays from 
both secondary images appear along a line in the sky since for each image the 
arrival directions move in opposite directions as the energy decreases below Eq. 
Then, this allows to measure the 'width' and 'length' (shortest and longest dis- 
tance) of the distribution of the events in the sky. In order to study the effect of the 
angular resolution, a random gaussian deflection with standard deviation a = 1° 
(which is the typical resolution of air shower arrays at ultra-high energies) is added 
to the position of each event. 

In Figure [6] we show the distribution of events within a certain width or length 
for sets of ten events from the one hundred simulated sources without and with 
angular resolution effects, performing in each case 10 4 different realizations of the 
arrival directions and energies of the events. As the typical width without exper- 
imental uncertainty is below 1°, the inclusion of the experimental angular uncer- 
tainty has a significant impact on this quantity, while in the case of the length, 
which for some sources is as large as 18°, the angular resolution effects have only 
a minor impact. From these distributions we obtain the characteristic angular di- 
mensions of the elongated excesses of events which would be seen in the sky due 
to the appearance of multiple images of a source. We use this to determine the 
angular window that should be adopted to perform the search of excesses due to 
multiple imaging. We will consider a maximum width of 2°, for which 92.5% 
of the simulated events are included (taking into account the angular resolution). 
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Concerning the maximum length, we set it at 8°, for which 92.8% of the sim- 
ulated events are contained. There is of course a compromise between a larger 
maximum length (or width), which will include the full length of the cosmic rays 
directions with energies between [0.9Eq, Eq] for a larger fraction of simulations, 
and the increase of the number of events expected from the background. We note 
that the longer structures are limited to certain regions of the northern folds near 
the galactic center. On the other hand, restricting the search to structures of shorter 
length would just imply loosing some events with energies close to 0.9 Eq and typ- 
ically belonging to the least amplified secondary image. We note that 86.5% of all 
the simulated events are located within angular windows of 8° length and 2° width. 

Having characterized the typical signal of clustering due to the appearance of 
multiple images of point-like sources, we construct an algorithm to search for this 
kind of features in cosmic ray data. The idea is to scan the sky counting the num- 
ber of events with energies in a range [0.9-E, E], within an angular window of 
2° x 8° arbitrarily oriented and compare the result with the expected number for an 
isotropic distribution of events, as a function of E. The significance of the excess 
is calculated with the Li & Ma method ETTl . For the simulations done in this pa- 
per, we will consider a full-sky observatory with equal exposure in all directions. 
Real observatories have an exposure which is not uniform on the celestial sphere 
and this fact can be easily taken into account and should not affect our results sig- 
nificantly. Due to the scan performed both in the location and orientation of the 
angular window as well as on the energy, in order to compute the chance probabil- 
ity of finding an excess with a given Li & Ma significance we just apply the same 
algorithm to simulations of isotropically distributed events to find what fraction of 
the isotropic simulations shows more clustering than the data. 



4 Application to simulations 

In this section we apply the method previously discussed to the simulations of 
sources with multiple images. To begin with, we analyze how many events an ex- 
cess would need to have in order that the probability for it to happen by chance 
be smaller than 10~ 3 . This is determined by analyzing simulations of isotropically 
distributed background events superimposed with events simulated from each of 
the sources with secondary images discussed in the previous section. The energy 
of the background events is selected following a broken power-law spectrum E^ 1 , 
for which we adopt the values determined by the Pierre Auger Observatory 11221 . 
7 = 2.59 for 4 EeV < E < 28.8 EeV and 7 = 4.3 for E > 28.8 EeV. For each 
of the simulated sources, the secondary images appear at a different rigidity Eq/Z 
and due to the falling cosmic ray spectrum, the number of background events de- 
pends on energy (and on the total number N to t of cosmic rays measured above 
the threshold). Therefore, the number of events from a source required to have a 
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Figure 7: Number of events with energies between [0.9 Eq, Eq] that are necessary 
to have an excess with isotropic chance probability smaller than 10~ 3 , as a function 
of the energy Eq at which the multiple images appear (a) and the corresponding 
fraction x of the cosmic ray flux in that energy range (b). 



detectable excess depends on Eq/Z and N to t- We present the results assuming the 
data set consists of Ntot = 10 4 measured cosmic rays above 9 EeV, comparable 
to the present data set of the largest existing air shower array, the Pierre Auger 
Observatory. The number of events with energy in the range [Q.QEq,Eq\ an d in 
an angular window of 2° x 8° coming from each of the one hundred simulated 
sources that are needed in order that the probability that an equal or larger excess 
appears by chance be 10~ 3 is plotted in Figure [7(a)] as a function of the energy of 
the caustic assuming proton composition of the source {Z = 1). The number of 
required events increases as the energy decreases (from 7 for Eq « 20 EeV to 1 1 
for Eq rj 10 EeV), which is due to the fact that the background signal grows at 
lower energies. The slight increase at energies above 29 EeV is due to the change 
in the slope of the broken power law spectrum. We note that if the total num- 
ber of ultra-high energy cosmic rays is increased, the number of events required to 
identify a secondary image with the specified significance will scale as \J iV toi /10 4 . 

Measuring such an excess would allow to identify a direction and an energy at 
which a caustic crosses the position of a source. This information would be valu- 
able in distinguishing between different galactic magnetic field models which lead 
to different predictions for the location of caustics. 

The possibility to detect the signal of secondary images depends on whether 
there is a strong enough source of cosmic rays in a region crossed by a caustic so 
as to produce the required number of events. This is directly related to the density 
of cosmic ray sources in our neighborhood since a larger source density implies 
fewer expected events from each source, what would make less probable the detec- 
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tion of this effect. 



An estimate of the maximum source density for which it should be possible 
to detect the events of the secondary images can be obtained by assuming equal 
intrinsic luminosity sources. We first determine for each simulated source the frac- 
tion x{Eq) of the cosmic ray flux in the energy range [0.9 Eq, Eq] that corre- 
sponds to the events that are necessary to detect the caustic with chance probability 
smaller than 10~ 3 (plotted in Figure |7(b)| ). This fraction rises from 4 x 1CT 3 at 
10 EeV to around 10~ 2 at 30 EeV. The flux from the secondary images in the 
energy range [0.9 Eq, Eq] is largely amplified with respect to the one in the ab- 
sence of lensing effects for most of the sources. Therefore, the required absolute 
luminosity of the source would be reduced by a factor l/// sec , with [i sec the ratio 
of the sum of the fluxes from the two secondary images in the energy range [0.9 
Eq, Eq] to the flux of the source that would be received in the absence of mag- 
nification effects in the same energy range (shown in Figure |8(a)| for each of the 
one hundred simulated sources as a function of their caustic energy). The fraction 
x^Eq) = x(Eq) I ' n sec {Eo) that would correspond to the required contribution of 
the source without lensing effects is plotted in Figure |8(b)[ For the sources in the 
regions of the sky where the caustic magnification is more important, this fraction 
turns out to be smaller than 10~ 3 . The prospects to detect the presence of caustics 
are best in these regions. From Figure |8(b")] it can be seen that x^(Eq) is smaller 
than 10~ 3 for 41 of the 100 simulated sources. This corresponds roughly to a frac- 
tion F = 1/10 of randomly chosen directions in the sky. The closest source in one 
of these regions is hence the best candidate to lead to a caustic detection. 

Assuming an homogeneous distribution of sources with mean density n and 
intrinsic luminosity L, we can estimate the mean flux expected from the closest 
source (/i) as 

= r 4^ ( eX P(-r F - 3 )) (™vrr 2 dr) = (Fnf* , 

(2) 

where the first parenthesis in the integral represents the probability of having no 
source within a distance r in that sector of the sky and the second one represents 
the probability of having one source between r and r + dr. 

On the other hand, the mean total flux from all sources in the sky located at a 
distance smaller than R m ax from Earth is given by 

{f tot) = LnR max . (3) 

Then, we can estimate the fraction of the mean flux from the closest source in a 
fraction F of the sky to the total flux as 

x (A) rq/3) f 2 / 3 

" {ftot) 62/3^1/3 n l/ 3jRmax - ^ 



12 



|-iaft ,■ T ,-, .-. - 

10 15 20 25 30 35 40 



164-1 1 1 1 1 — . — i 1 1 — . — i 1 — 

10 15 20 25 30 35 40 



EJEeV] EJEeV] 
(a) (b) 

Figure 8: The ratio of the sum of the fluxes from the two secondary images in 
the energy range [0.9 Eq, Eq] to the flux of the source that would be received in 
absence of magnification effects in the same energy range (/x gec ) (a) and the fraction 
Xn(Eo) that would correspond to the required contribution of the source without 
lensing effects (b) as a function of the critical energy Eq at which the multiple 
images appear. 



This allows to estimate an upper limit on n given the fraction of the flux x^ 
necessary to detect the signal from the secondary images and the fraction F of the 
sky directions. For example, for x M = 10~ 3 and F = 0.1 and considering that cos- 
mic rays of energy E ~ 10 EeV can arrive from R ma x — 1 Gpc, the upper limit 
on the source density for which it is expected to be possible to detect an excess 
due to secondary images of a source is n ~ 2 x 10~ 3 Mpc~ 3 . This large value 
is compatible with the lower bounds that could be deduced from the searches of 
clustering at the highest energies. 

This is a rough estimate of the mean density value with the simplifying as- 
sumption of equal intensity sources with an homogeneous distribution, so large 
fluctuations of this value can be expected. However, the order of magnitude of the 
obtained result implies that there are reasonable chances to detect this effect in the 
near future, at least if there is a nearby source of light ultra-high energy cosmic 
rays. For a source of heavy nuclei, the caustics will appear at energies Z times 
larger and hence the observation of these effects will be much harder. 



5 Conclusions 

Magnetic lensing effects on ultra-high energy cosmic rays produce interesting phe- 
nomena which could be measured with the increasing statistics of present and fu- 
ture air shower experiments. A remarkable aspect is the appearance of multiple 
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images of a source, which could be detectable due to the magnification of the flux 
near the critical energy at which the secondary images appear. 

In this work, the typical angular distribution of secondary images was charac- 
terized with simulations of sources that present multiple images in a BSS-S model 
for the galactic magnetic field. Consequently, an algorithm to search for clustering 
of events of similar energy in angular stripes of 8° x 2° in cosmic rays data was im- 
plemented. The significance of the excess is calculated with the Li & Ma method 
and the probability of finding such a significance is determined from isotropic sim- 
ulations. 

We determined the minimum number of events that are required to detect an 
excess with probability smaller than 10~ 3 of occurring by chance in an isotropic 
distribution of arrival directions. We estimated also an upper limit to the source 
density which is needed to be able to detect such a number of events from a source. 
Our results indicate that such an excess should be detectable in the near future with 
the current ultra-high energy cosmic rays observatories if the source density is not 
too large and if there are nearby sources of light ultra-high energy cosmic rays. 
Detecting an excess due to the appearance of multiple images of a source would 
also allow to obtain information regarding the location of caustics which would be 
valuable in disentangling between different models for the galactic magnetic field. 
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